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Abstract 

We present a detailed study of the optical absorption spectra of DNA bases and base pairs, 
carried out by means of time dependent density functional theory. The spectra for the isolated 
bases are compared to available theoretical and experimental data and used to assess the accuracy 
of the method and the quality of the exchange-correlation functional: Our approach turns out to 
be a reliable tool to describe the response of the nucleobases. Furthermore, we analyze in detail 
the impact of hydrogen bonding and vr-stacking in the calculated spectra for both Watson-Crick 
base pairs and Watson-Crick stacked assemblies. We show that the reduction of the UV absorption 
intensity (hypochromicity) for light polarized along the base-pair plane depends strongly on the 
type of interaction. For light polarized perpendicular to the basal plane, the hypochromicity effect 
is reduced, but another characteristic is found, namely a blue shift of the optical spectrum of the 
base-assembly compared to that of the isolated bases. The use of optical tools as fingerprints for 
the characterization of the structure (and type of interaction) is extensively discussed. 
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1. INTRODUCTION 



Optical absorption and circular dicliroism are widely used experimental tools to charac- 
terize the structural and dynamical properties of biomolecular systems. The particular merit 
of optical characterization tools is that they can discriminate between intrinsic molecular 
properties and solvent-induced effects. This is especially important for DNA and DNA-based 
compounds.i^ Traditionally, DNA molecules have always retained a special place in scien- 
tific investigation, for biological/medical issues. Lately, DNA is also attracting interest for 
several potential applications in the field of nanotechnology, due to its stability (in solution), 
to its one-dimensional character, and to the regular vr-stacking, along with the unique proper- 
ties of self-assembling and recognition . '^'^'^i'^'^ In relation to nanotechnology exploitation, the 
determination and interpretation of the electronic properties of nucleobases and of DNA he- 
lical arrangements are an extremely valuable foreword, and notable multidisciplinary efforts 
are currently devoted to such goals: we refer the reader to recent reviews about electronic 
structure calculations and possible charge motion behaviors.—"^ Furthermore, the knowledge 
of the electronic properties, excited state lifetimes, and ultra-violet (UV) absorption spectra 
is of paramount importance for our understanding of, e.g., the crucial phenomenon of UV 
radiation-induced DNA damage . From this brief preamble, it is clear that the full 
characterization of the optical properties of DNA molecules and DNA-based complexes is 
of great interest. In order to relate the optical properties of the nucleic acids to their struc- 
ture, spatial conformation, and type of intra-molecular interactions, a valuable preliminary 
step is to gain insight into the excited-state properties of their building blocks, namely the 
monomeric bases, and to understand the role of hydrogen-bonding and stacking when these 
monomeric units form complex assemblies. In their natural environment, the DNA bases 
are paired via hydrogen-bonds in the Watson-Crick scheme,— and are covalently bonded to 
the sugar-phosphate backbone. The hydrogen-bonded base pairs interact with each other in 
the typical helical arrangement by inter-plane van der Waals forces.— To disentangle how 
the different interactions control the DNA dynamics upon light absorption, it is important 
to infer how the spectrum of a given isolated nucleobase is modified by mutual interactions 
in the different spatial conformations of DNA-assemblies. This is the goal of the present 
work: To provide a systematic study of the stacking and H-bonding interaction effects in the 
optical spectra of molecular complexes formed with isolated nucleobases. We undertake this 
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task by means of time-dependent-density functional theory (TDDFT) calculations.— li^ii^^ 
TDDFT is gaining increasing popularity as an efficient tool to calculate electronic exci- 
tations in finite systems, thanks to its simplicity and moderate computational cost. Since 
the very beginning, TDDFT has provided good results for the optical response of a large 
set of molecular systems,-^ including some biomolecules (e.g., green fiuorescent protein chro- 
mophores and their mutants , chlorophylls,— fiavins^i and nucleic acid bases , among 
many others). It is still perhaps premature to discuss the general level of accuracy of TDDFT 
when applied to biomolecules, especially when dealing with van der Waals complexes."^! Fur- 
thermore, as the field is quite recent, we can still expect rapid methodological developments, 
in particular toward the derivation of better exchange- correlation functionals. However, the 
existing results are already promising, and pave the way for a broader application of TDDFT 
in biochemistry. 

Before studying complex helical DNA-based biopolymers, extensive tests are required to 
prove the performance and predictive power of TDDFT for DNA-based systems. One of 
the aims of the present work is to provide results from first-principles calculations of simple 
DNA-based assemblies that can be used as reference for future developments and studies. 
This reference set includes the isolated DNA bases, hydrogen bonded Watson-Crick pairs, 
a stacked guanine-cytosine dimer, and a quartet formed by two stacked guanine-cytosine 
pairs. 

For the isolated DNA bases there are plenty of experimental and theoretical results con- 
cerning their optical response.— We do not carry out a full systematic analysis of all the 
published data for the isolated DNA-bases, but only those more pertaining to our investiga- 
tion. The computational tools applied to nucleobases range from single excitation configura- 
tion interaction (CIS), to complete active space 2nd-order perturbation theory (CASPT2), 
and include other approaches such as multireference perturbation configuration interaction 
(CIPSI). A review of the performance of such methods on the nucleobases,— as well as rele- 
vant experimental data can be found in Refsl^l^l^^- 

Despite the large number of works on the excited states of isolated DNA bases, a lim- 
ited number of studies is devoted to base pairs and base assemblies due to their complex- 
ity. Shukla and Leszczynski studied adenine-uracil,— adenine-thymine (AT) and guanine- 
cytosine (GC) pairs in the Watson-Crick configuratioii^^ at the CIS level; Sobolewski and 
Domcke^ studied the low lying energy part of the spectrum of GC base pairs with the more 
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sophisticated CASPT2 technique; Wesolowski^ used an embedding method to study the 
lowest excited state of the GC and AT pairs. Very recently, Tsolakidis and Kaxiraa^^ com- 
puted the whole absorption spectra of the GC and AT pairs in different tautomeric forms 
of the bases in the TDDFT framework. 

Moreover, very few studies exist on the excited-state properties of bases in a stacked 
configuration. Jean and Hall studied the fluorescent properties of dimers of 2-aminopurine 
stacked with DNA bases^- in different forms and stacked trimers containing 2-aminopurine,— 
showing the relevance of the stacked geometries in the character of the excited state tran- 
sitions. To the best of our knowledge, the present results are the first ab-initio calculations 
dealing with stacked natural pairs. Yet, we remark that the present results are for free 
standing nucleobase complexes, i.e., not including solvation effects. These effects are known 
to be more important for nn* than tttt* type transitions . 

The paper is organized as follows: In Sect.|21we provide the details of the real-time real- 
space method^ used to compute the optical properties. In Sect. El we present our results 
for the isolated gas-phase nucleobases, i.e., guanine, cytosine, adenine, thymine, and uracil, 
(G, C, A, T, and U, respectively). In order to understand the role of hydrogen-bonding 
in shifting and modifying the spectral features, we show, in Sect. |31 results for H-bonded 
Watson-Crick GC and AT base pairs (labeled GCh and ATh, respectively). In Sect. El we 
present results for a GC stacked dimer (labeled GCs) that mimics the arrangement between 
C and G in two consecutive planes in the real DNA double helix. We also discuss the 
relative roles of 7r-stacking and hydrogen bonding in the optical absorption of a stacked 
quartet made of two adjacent Watson-Crick GC pairs as in A-DNA [labeled d(GC)], that 
combines hydrogen-bonding and vr-stacking. In Sect. El we summarize the results of the 
present work and provide some perspectives for future studies. 

2. COMPUTATIONAL FRAMEWORK 
2.1. Optical spectra by TDDFT 

The optical absorption spectra of DNA bases and base pairs were computed within a real- 
space real-time version of TDDFT, as implemented in the code octopus.— This method does 
not rely on perturbation theory and is competitive with implementations in the frequency 
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domain The theoretical background and the computational details of our scheme 
are extensively described elsewher o^^ ,42^46^47 . i^qj-q -^g j^g^ summarize the crucial aspects 
relevant to the present study of DNA complexes. 

The starting point is the calculation of the ground state electronic structure in the DFT 
framework, which is done within a pseudopotential approach (see below for details). In order 
to access the excited state properties, the ground state is then instantaneously perturbed 
with an electric field of magnitude along the three principal Cartesian directions (i.e., 
by applying the external potential v{r,t) = —kQXy5{t), where = x,y,z). The amplitude 
ko must be sufficiently small in order to keep the system dipolar response linear. In this 
impulsive approach, all the frequencies of the system are excited with the same weight. 
Next, the time dependent Kohn-Sham (KS) orbitals are propagated for a finite time, and 
the dynamical polarizability is obtained as: 

1 ^.3. 



a^{uj) = -— d r Xu6n{r,uj) , (1) 

Kq J 

where Sn{r,uj) is the Fourier transform of the relative density n{r,t) — no{r), and no(r) 
stands for the ground state density. The photo-absorption cross-section, the quantity usually 
measured in experiments, is proportional to the imaginary part of the polarizability averaged 
over the three space directions: 

^M = ^^EM^), (2) 

where c is the velocity of light. Another widely used quantity is the dipolar strength function 
Siy^Lj)^^^^ which is connected to a^{uj) by 

Su{lj) = ^^'^ai,{uj) . (3) 

With this definition, the Thomas- Reiche-Kuhn sum rule integrates to the number of electrons 
in the system. 6*^(0;) is the quantity plotted in our figures. 

A significant advantage of this real-time approach is the fact that only occupied states 
are needed, thus avoiding the calculation of the unoccupied states that instead enter the 
traditional orbital (occupied-empty) representation of the linear response equation.^^ The 
complete set of empty orbitals required in the latter approach is fully accounted for by 
the time propagation. In addition, it is important to note that in the time domain only 
the approximation to the exchange-correlation potential v^c is required, whereas in con- 
ventional frequency-domain linear response TDDFT formulations also the f^c kernel is 
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Figure 1. Dipole strength function of cytosine along the axis perpendicular to the base 
plane, calculated through the real-time propagation (blue solid curve) and in the frequency 
domain^ (green dashed curve). An artificial broadening of O.lSeV was adopted as an eye 
help to compare time and frequency domain spectra. Note that we can distinguish the nvr* 
character in an indisputable way only for the lowest energy peak. At higher energy the nn* 
character is less pure, due to the mixing of a large number of transitions. This fact is 
illustrated here for cytosine, but we observed it for all the investigated nucleobases. 

needed . ^^i^^i^'^'i'^^i'^^i'^^i"^^ For v^c we employed the local density approximation (LDA) with 
the Perdew-Zunger parametrization.^^ In the case of the isolated guanine base, we also per- 
formed the time propagation with the generalized gradient approximation (GGA) PW91^ 
fxc functional. As expected, no substantial differences were detected between the two pa- 
rameterizations. Therefore, and as the LDA functional is numerically the most stable, we 
chose to use this functional in the time propagation for all the other investigated systems.— 

Concerning technical details: We used a real-space grid made of overlapping spheres 
with a radius of 4 A centered around each nucleus. All quantities were discretized in a 
uniform grid with spacing 0.23 A. The time-step for the time evolution was 0.0066 fs (which 
ensures stability in the time-dependent propagation of the KS wave functions) and the 
total propagation time was at least 20 fs. Such simulation parameters ensure well converged 
absorption spectra up to about 8eV. The energy resolution, dictated by the total simulation 
length by AE ~ h/T, is better than 30meV. This sets the lower limit for the linewidth 
of the calculated spectra. The electron-ion interaction (both in the time-dependent and 
time-independent DFT calculations) was modeled by norm-conserving pseudopotentials.-^ 

To compare with experiments where the molecules are not aligned with respect to the 
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applied field, we need to average the response along the three Cartesian axes [cf. Eq. Q]. 
This result can be directly compared to other available data and used to discriminate among 
different nucleobase conformations. However, the analysis can be pushed even further: in 
the case of planar molecules (isolated bases, GCh, ATh) one can easily distinguish between 
in-plane and out-of-plane transitions. While the former are prevalently tctt* and have large 
oscillator strengths, the latter are much weaker (by one or two orders of magnitude) and 
are usually hidden under the stronger signals. In order to unveil the character of the weak 
transitions, the time-domain analysis is not sufficient: Thus, we also performed linear re- 
sponse calculations in the frequency domain*^^ for all planar structures (isolated bases, GCh, 
and ATh pairs). In this formalism, several pairs of occupied/unoccupied orbitals participate 
in each excitation with different weights: it is then possible to qualitatively assign to each 
excitation the predominant character of the contributing transitions. In Fig. ^we show an 
example of the transition character of spectral features, derived from our frequency analysis 
of cytosine. The dipole strength function of cytosine along the direction perpendicular to the 
base plane, obtained by both the time-domain and frequency-domain techniques, is plotted 
against energy. The agreement encountered between the two techniques turns out to be 
good in the accessible energy range. For all other planar systems, we present the analysis 
of the transition character obtained in a similar way (namely, by comparing time-domain 
and frequency-domain features), but without explicitly showing the figures for the spectra 
obtained by frequency-domain TDDFT. 

We note that in the frequency-domain the transitions are delta functions. In order to 
make a fair comparison with the time-domain propagation, we have artificially broadened 
both spectra in Figure 1 using a lorentzian function of 0.15 eV width. Note that for the rest 
of the spectra shown in the paper obtained with the time-propagation scheme, we did not 
add any additional broadening and the resolution of 30meV is fixed by the duration of the 
time evolution. Any additional broadening comes from the coupling of the excited state to 
other excitations (e.g. Landau damping). 
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2.2. Ground-state equilibrium configuration by DFT 

1. Starting atomic configurations 

The starting atomic coordinates of the isolated DNA bases and of the ATh pair were 
derived from standard structural parameters for DNA,i^ whereas for the systems GCh, 
GCs, and d(GC) we extracted the simulated fragments from a high-resolution X-ray crystal 
structure of A-DNA,^^ after a classical force-field structural adjustment. The crystal is 
formed by short helices of double stranded DNA with the sequence d(AGGGGCCCCT), a 
model for poly(dG)-poly(dC) DNA. The H coordinates (absent in the X-ray PDB file) were 
assigned using the HBUILD command^ of the CHARMM package.^^ The CHARMM force 
field was also used for a raw relaxation of the entire d(AGGGGCCCCT) structure. From 
this grossly relaxed polymer we extracted the H-bonded and stacked GC fragments, that 
have been subsequently relaxed by quantum simulations. The simulations allow us to have a 
complete set of base configurations that could be used to perform configurational sampling 
and extract structural broadening effects, the analysis of this issue is left for future work. 

All our simulations were done for systems in the gas phase: Note that the effect of a sol- 
vent may be crucial when comparing the computational results with experimental data.>22i^ 
Furthermore, the sugar-phosphate backbone was neglected. This restriction should not be 
relevant in the energy range we are interested in (3-7.5 eV), as the sugar and the phosphate 
contribution to the absorption spectra only starts to be important above 7.0 eV. 

2. Technical details of the DFT runs 

The DFT structural optimization was carried out using the B3LYP/6-31++G(d,p) 
exchange-correlation functional for the isolated bases and GC assemblies, and the PW91 
functional^^ for the ATh pair. This latter functional was widely tested in guanine-rich 
systems,— i^i*^ and yields a similar accuracy for ATh pairs and stacked pairs. 

It is well known that LDA may be grossly wrong in the description of H-bonds: that is why 
different gradient-corrected functionals were used to relax the geometries. However, we also 
know that, once the geometry is given, the the LDA and GGA furnish very similar results for 
the excitations, therefore we stick to the LDA functional for simplicity. We are confident this 
is a reasonable approximation. Indeed, in the test mentioned above for guanine, different 
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Exp. TDDFT(^) TDDFTC^) CASPT2('=) TDDFT^'^) 
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7.61 
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TABLE I: Vertical excitations energies (eV) calculated for cytosine, compared with 
averaged experimental values and selected computational results. The columns correspond 



to: (a) this work; (b) Ref, 



25; fc) Ref 



SOt (d) Refl26l The experimental absorption and CD 



iorpt 



values, along with the original references, are collected in Reft 
functionals do not yield significant differences as far as the optical excitations are concerned. 



3. ISOLATED GAS-PHASE NUCLEOBASES 



In this section we present the calculated absorption spectra for the five isolated nucle- 
obases. Both the purines and the pyrimidines exist in nature in different tautomeric forms. 
We limited our calculations to the 9H keto form of guanine, and to the 9H amino form of 
adenine. The latter tautomer is the one present in DNA and RNA polymers, and is therefore 
the most relevant regarding the adenine properties in nucleic acids. 



3.1. Pyrimidines 



1. Cytosine 



The computed total and perpendicular photo-absorption cross sections of cytosine are 
shown in Fig. El The energies of the spectral peaks are reported in Table d together with 
other selected sets of data that represent the state-of-the-art performance of experimental 
and computational techniques on cytosine. 
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Figure 2. Photo- absorption cross section of isolated gas-phase pyrimidine nucleobases. 

The solid blue (dashed yellow) line is the signal averaged along the three real-space axes 
(projected onto the axis perpendicular to the base plane). Insets: HOMO (left) and 
LUMO (right) Kohn-Sham wave functions. The cyan (magenta) isosurfaces represent 
positive (negative) charge values. Different atoms are indicated with different colors: 
carbon (green), nitrogen (blue), oxygen (red), hydrogen (white). 



ance at Table IJ reveals a rather 
25l . and a poorer matching with 



For what concerns the transition energies, an overall gl 
good agreement of our results with those obtained in Ref 
the outcome of B3LYP-TDDFT based on GAUSSIAN98,^ and of CASPT2. None of the 
computational data is in perfect accord with experiment. Therefore, a selection of the 
best theoretical approach on the basis of the comparison between theory and experiment is 
hindered. Our agreement with experimental values is much more satisfactory in the high- 
energy than in the low-energy regime. 
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Figure 3. The same as Fig. |21for adenine and guanine. 



If we turn our attention to the relative oscillator strengths between the various peaks, 
then the agreement with the data of Ref|2^is not so perfect: we encounter differences in the 
distribution of oscillation strengths in the high-energy range of the spectrum. We will see 
later that these high energy discrepancies also appear in the spectra of the other nucleobases 



25, and absent in 



and can be traced to the limitations of the localized basis set used in Ref 
the present work due to the use of a real-space grid. 

Combining the time-domain spectra with the frequency- domain analysis, we can describe 
the spectral peaks in terms of transitions between electron states. Fig. Q shows an nir* 
transition at 4.3 eV in the perpendicular spectrum, not clearly discernible in Fig. El because 
of the low intensity of the perpendicular signal. This means that the lowest energy peak 
at 4.18eV in the averaged spectrum (Table H)), absent in the perpendicular polarization, 
must have a vrvr* character. The fact that the first optical excitation in cytosine is vrvr*- 
like was predicted in all recent computational studies. The spacing between the lowest vrvr* 
excitation and the mr* transition computed by us is 0.12 eV: this value is in good agreement 
with the results of BSLYP-TDDFT^-S^ and CIS,-^ whereas CASPT2 calculations^^ predict 



a larger spacing between the two excitations. A smaller spacing was found in Ref 
B3LYP-TDDFT. 



65 with 
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Exp. TDDFT(^) CASPT2('^) TDDFT^'^) 



4.6-4.9 


4.69 


4.54 (nvr*) 


4.66 (nvr*) 


4.9-5.2^ 


5.25 


5.00 


5.17 


5.8-6.1 


5.92 


5.82 


5.89 


6.3-6.6 


6.28 


6.46 


6.42 


6.8-7.0 


6.98 


7.00 


6.81 




7.57 







TABLE II: Vertical excitations energies (eV) calculated for uracil, compared with 
average experimental values and selected calculations. The columns correspond to: (a) this 
work; (b) RefH (c) RefQ; (d) this band was obtained in CD spectra and it was argued 
that it exhibits nvr* character. The experimental absorption and CD values, along with the 



original references, are collected in Refj28l: the second experimental band that is taken 



67 



68. 



from Refs 

2. Uracil and thymine 

The computed total and perpendicular photoabsorption cross sections of uracil and 
thymine are shown in Fig. El The energies of the spectral peaks are reported in Tables ITTl 
and mil The spectra of the two bases U and T are very similar, both in the number of peaks 
and spectral energies. A fairly good agreement is met with averaged experimental data, as 
well as with TDDFT^S^ and CASPTS^^ computational data. 

All previous calculations, either in the TDDFT framework or with quantum chemistry 
methods, predict the lowest transition to have an mr* character in vacuo. We agree with this 
assignment for uracil (the perpendicular signal is not visible at 4.69 eV in Fig. |21due to the 
weak intensity). In the case of thymine, we do not find any appreciable signal perpendicular 
to the base plane contributing to the peak at 4.54 eV. However, this is consistent with 
TDDFT— and CASPT2^ calculations, as they indeed report nir* transitions with extremely 
small oscillation strengths (le-4-le-6), unresolved by us. 

The character of the band at 4.9-5.2 eV for uracil and at 5.0-5.1 eV for thymine is not yet 
clear: the absorption spectra of uracil compounds^ and circular dichroism experimentsSi in- 
dicate that it is due to a vrvr* transition. However, such vrvr* character has not been confirmed 
by magnetic circular dichroism^^ and polarized absorption spectra experiments.— i^i^^^iZ^ 
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Exp. TDDFT(^) TDDFTC^) CASPT2('=) TDDFT^'^) 



4.5-4.7 


4.54 


4.45 


4.39(n7r*) 


4.69 


5.0-5.1 


5.30 


5.24 


4.88 


4.96 


5.8-6.0 


5.74 


5.68 


5.88 


5.95 


6.3-6.6 


6.36 


6.38 


6.10 


6.19 




6.51 






6.50 


7.0 


6.88 


6.86 


7.13 


6.86 






7.07 








7.21 


7.28 








7.59 


7.52 







TABLE III: Vertical excitations energies (eV) calculated for thymine, compared with 
average experimental values and selected calculations. The columns correspond to: (a) this 
work; (b) RefQ (c) RefQ (d) Refl The experimental absorption and CD values, along 

with the original references, are collected in R,ef ]28l: the experimental second band that is 

taken from Refs|i 

Lorentzon and coworkers, by performing calculations at the CASPT2 level,— suggest that 
this band has a nvr* origin: they arrived at this conclusion by correcting their computed 
values with a 0.5 eV blue-shift due to the solvent. A similar conclusion was reported for 
uracil by Shukla and coworkers,^^ who performed CIS calculations taking the solvent into 
account with a polarizable continuum model. However, more recently, Shukla and Leszczyn- 
ski using B3LYP-TDDFTi2& reported only a tttt* transition at 5.17eV for uracil and at 4.96 
for thymine in the energy range under consideration. In our work, we do not detect any 
absorption in the direction perpendicular to the base-plane in this energy range, whereas 
peaks at 5.24 eV for T and at 5.25 eV for U are found, with considerable oscillator strength 
induced by light polarized in the plane of the molecule. Therefore, our results indicate that 
the second band most likely has a tttt* character. 
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Exp. TDDFT TDDFTC^) CASPT2('=) TDDFT(^) 



4.59 


4.51 


4.51 


5.13 


4.94 


4.8-4.9 


4.88 


4.95 


5.20 


5.21 


5.38 


5.49 


5.58 






5.7-6.1 


5.72 


5.79 


6.24 


5.93 


6.26 


6.19 


6.28 


6.72 


6.12 




6.49 


6.63 




6.16 


6.81 


6.76 


6.92 


6.99 






7.09 


7.47 






7.73 


7.39 


7.81 


7.57 





TABLE IV: Vertical excitations energies (eV) calculated for adenine, compared with 
averaged experimental values and selected calculations. The columns correspond to: 
(a) this work; (b) Ref 2^; (c) RefQ; (d) Ref|3 The experimental absorption and CD 
values, along with the original references, are collected in R,efl29[ 



3.2. Purines 



1. Adenine 



The computed total and perpendicular photo-absorption cross sections of 9H-adenine 
are shown in Fig. El The spectral peaks are summarized in Table IIVI and compared with 
other available theoretical and experimental data. An excellent agreement with other LDA- 
TDDFT computed valuesSS= is observed, regarding both the peak energies and the relative 
oscillation strengths. Slight discrepancies occur only in the high energy range: these are 
most likely due to differences in the technical details adopted in the calculations (e.g., the 
basis sets). 

Looking at the spectrum in the direction perpendicular to the base plane, we can dis- 
tinguish four prevalently mi* transitions in the energy range 4.10-4.79 eV, that are hidden 
in Fig. 121 because of their tiny oscillator strengths. We find that the lowest excitation has 
mainly nvr* character and is very close in energy to the first tttt* transition (forming the 
peak at 4.51 eV). This result is in agreement with resonant two-photon ionization and laser 
induced fluorescence spectroscopy of jet cooled adenine,'^^ provided that the adiabatic tran- 
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sition energies follow the same trends as the calculated vertical ones. Our findings also agree 
qualitatively with frequency-domain B3LYP-TDDFT— and CIPSSLZ^ calculations, whereas 
CASPT2^^ and CIS^>^S^ yield the reverse order. Given the small energy difference between 
the two transitions, we can conclude that the overall qualitative agreement is satisfactory. 

Proceeding to higher energies, we find the second vrvr* peak 0.37 eV after the first one. 
Experiments indicat o^^i^^i^^ that the low energy portion of the adenine photo-absorption 
spectrum consists of two closely spaced tttt* transitions. Holmen and coworkers showed by 
linear dichroism that those two peaks are separated by 0.26 eV.— This value is in reasonable 
agreement with our results and with other TDDFT calculations,-^^*^ while different tech- 
niques predict a closer energy spacing. On the other hand, measured spectra indicate that 
the lowest energy transition has a smaller oscillator strength than the subsequent one: such 
an evidence is consistent with the CASPT2 and CIPSI results, whereas TDDFT simulations 
yield an inverse ordering of oscillator strengths in the first two peaks^^ (see Fig. E)). 

Concerning other rnr* transitions, there is experimental evidence by circula r^^i^° and 
lineal^ dichroism of a signal around 5.4 eV. Our calculation is also able to reveal a signal 
in this range. However, the agreement with the experiment for the nvr* transitions should 
be considered with caution, because of their high sensibility to the effect of the solvent 
(neglected in our calculations) and exchange-correlation functional. 



2. Guanine 

To conclude the presentation of the results for the isolated DNA bases, we show in Fig. El 
the computed total and perpendicular photo-absorption cross sections of 9H-guanine. The 
spectral peaks are summarized in Table |3 and compared with other available data. The 
excitation energies are in good agreement with the averaged experimental data and with 
the results of Refl^- Small differences are encountered in the high energy region of the 
spectrum in the relative oscillator strengths of the peaks at 6.22 eV and 6.67 eV, while the 
low energy region is in perfect agreement. Again, we attribute the discrepancies at higher 
energies to the different basis sets used in Ref|3 and in our work. In the region from 5 to 
6 eV we also find very weak peaks that are not observed experimentally. 

Regarding the out-of-plane spectrum, we find one peak at 4.47 eV that has a purely nir* 
character. Thus, the first mr* transition has a higher energy than the first tttt* transition 
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Exp. TDDFT TDDFT'^ CASPT2'^ TDDFT"^ 



4.4-4.5 


4.40 


4.46 


4.76 


4.85 


4.9-5.0 


4.66 


4.71 


5.09 


5.11 




5.01 


5.04 








5.28 








5.7-5.8 


5.76 


5.64 


5.96 


5.59 


6.1-6.3 


6.22 


6.23 


6.55 


5.83 


6.6-6.7 


6.67 


6.53 


6.65 










6.66 








6.82 


6.77 






7.04 


6.93 








7.58 


7.26 







TABLE V: Vertical excitations energies (eV) calculated for guanine, compared with 
averaged experimental values and selected calculations. The columns correspond to: 
(a) this work; (b) Refj^ (c) Refl3 (d) Ref 26. The experimental absorption and CD 



values, along with the original references, are collected in Ref 
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found at 4.40 eV (see Table \V\ . This attribution is consistent with previous CASPT2,— 
CIS,.2a and BSLYP-TDDFT^^ calculations. 



3.3. General comments on the results for the isolated gas-phase nucleobases 

The results of our time- dependent calculations^ of the excitation spectra of DNA bases 
and uracil in the gas phase show a satisfactory agreement with the experimental data and 
with other computational approaches, especially for the most intense peak above 6eV. In 
particular, good agreement is found with recent LDA-TDDFT calculations using a localized- 
basis set.'^^ Minor differences are found regarding the distribution of the oscillator strengths 
in the high energy range of the spectra. Such discrepancies can be ascribed to the different 
basis sets: a uniform grid in real space in this work, and a set of localized atomic orbitals 
in Ref H 

Regarding the more problematic nvr* transitions, we find that the weak perpendicular 
signal is sometimes hidden below the in-plane signal. In the cases where the nvr* transi- 
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tions can be detected, we find the correct relative ordering between nir* and tttt* peaks, as 
compared with other TDDFT calculations performed in the frequency domain. 

Another interesting trend is the relative ordering among different nucleobases for what 
concerns the first vrvr* transition. For this feature, supersonic jet experiments report C < 
G < ^.J4iS2i2i We indeed find that ordering in our first-principles calculations, i.e. the first 
excitation at 4.18eV for cytosine, 4.40 eV for guanine and at 4.52 eV for adenine. 

Last, we note that our calculations give the expected character of the Kohn-Sham orbitals 
of the different bases (see insets in Figs. |2]and|ni): the LUMO has vr character for all the 
nucleobases, whereas the HOMO is vr-like for the purines and cr-like for the pyrimidines, 
with a high charge density around the oxygen atoms in all cases. 

4. WATSON-CRICK PAIRS GCh AND ATh 

We now turn the discussion to the standard assembly of DNA bases, namely the Watson- 
Crick pairing. Figure ID shows the calculated TDDFT spectra of H-bonded GCh and ATh 
pairs (left) and the spectra resolved in the direction perpendicular to the basal plane (right). 
The sum of the photo-absorption cross sections from the isolated component nucleobases is 
also shown by the dashed lines (indicated by G+C and A+T), to reveal the effects induced 
by H-bonding in the Watson-Crick pairs. The peaks observed in the spectra are also listed in 
Table for a more detailed analysis. Figure El illustrates the isosurface plots of the HOMO 
and LUMO, which are the single-particle orbitals more important for the optical transitions 
in the energy range discussed here. In both pairs, the HOMO is localized at the purine and 
the LUMO at the pyrimidine. The HOMO-LUMO gaps are L93 eV and 3.07 eV for the 
GCh and ATh pairs, respectively. 

We focus our analysis on the most evident features induced by H-bonding base pairing, 
that are: (i) a small shift found for the lowest frequency peaks; (ii) the hypochromicity 
(intensity decrease) at high frequencies; (iii) an overall blue-shift of the spectrum for light 
polarized perpendicular to the pair. 

The relaxed structure of the GCh pair is practically planar. The shape of the total photo- 
absorption spectrum of the hydrogen-bonded GCh pair remains rather similar to the linear 
superposition of the isolated bases (cf. the top left panel of Fig. |3|). The first vrvr* peak at 
4.29 eV is composed of the first excitation of cytosine and the first excitation of guanine. We 
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Figure 4. Photo-absorption cross section of the GCh and ATh base pairs averaged (solid 
blue, left) and in the direction perpendicular to the base plane (solid red, right). The 
linear combinations of the spectra of the isolated purine and corresponding pyrimidine 
(G+C and A+T; dashed lines) are also shown for comparison. 

note a slight redshift of this peak with respect to the first peak in the red curve: this effect is 
the combination of a blue-shift by about 0.1 eV of the first cytosine peak and a comparable 
redshift of the first guanine peak. Such finding is in qualitative agreement with a CASPT2 
calculation by Sobolewski and Domcke^^ and by a calculation performed by Wesolowski^^ 
by embedding methods. 

For the second vrvr* transition at 4.72 eV the situation is reversed: The peak is blue-shifted 
with respect to the superposition of isolated G and C. This shift results from a blue-shift 
of the second guanine peak and a red-shift of the second cytosine peak. In the high-energy 
spectral range, we remark that the effect of the H-bonding pairing does not change the 
position of the brightest peak at 6.17eV, which has both G and C components. However, 
the intensity of this peak is depressed (hypochromism) by about 30% with respect to the 
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TDDFT(^) TDDFT(b) CASPT2('=) 


AT 

TDDFT TDDFTC^) 


2.44 


2.37 








3.14 


3.53 






3.22 


3.94 






4.05 


4.05 


4.29 


4.21 


4.35 




4.26 


4.72 


4.75 


4.67 


4.56 


4.40 






4.75 


4.98 


4.80 


5.38 


5.52 






5.20 


5.72 


5.69 






5.45 




5.84 




5.53 


5.60 




6.07 




5.87 


5.76 


6.17 


6.33 




6.20 


6.24 


6.46 


6.47 






6.48 


6.78 


6.68 




6.78 


6.73 




6.83 




7.80 


6.97 


7.07 


7.18 








7.53 


7.28 









TABLE VI: Vertical excitation energies (eV) for Watson-Crick GCh and ATh pairs, 
compared to other calculations. The columns correspond to: (a) this work; (b) RefQ 
(c) Refls^l (cf. also the results of Shukla and Leszczynski)^^; (d) Charge transfer transition. 

free monomers. Above this energy the spectrum results much more changed. 

Looking at the direction perpendicular to the plane of the GCh pah (top right panel), we 
observe an overall blue-shift of the spectrum. In particular, the first peak is shifted upward 
by 0.48 eV. This behavior reflects the well known uplifting of nir* excitation energies in 
hydrogen-bonding environments, already pointed out by other authors , ^^'^'^ and appears 
naturally from our first-principles simulation. 

For the ATh pair. Table IVII highlights an overall good agreement with a recent LDA- 
TDDFT calculation.^^ Small differences can be reasonably imputed to the different structure 
obtained upon relaxation. We find an almost planar geometry, in agreement with Shukla 
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and Leszczynski,~ while Tsolakidis and Kaxiras^^ find a propeller structure. 

As in the case of the GCh pair, the overall spectrum of ATh is quite similar to the linear 
combination of the spectra of the isolated A and T bases. A small blue-shift, by ~0.1eV, 
is detected in the low-energy range (roughly below 6 eV) with respect the superposition of 
isolated bases. In the high-energy part of the spectrum, we observe an enhancement of the 
strength with respect to isolated adenine and thymine for the peak located at 6.78 eV, and 
a reduction for the peak at 7.18 eV present in the isolated thymine spectrum. Hence, for 
the most intense peak at 6.78 eV, the effect of H-bonding is contrary to the hypochromicity 
described above for the GCh pair. 

In the out-of-plane spectrum (bottom right panel of Fig. 0}, as for the GCh pair, we 
observe that the effect of hydrogen bonding is an overall blue-shift. In particular, the first 
peak (that has an nn* character) is shifted by 0.4 eV; the second and fourth peaks are 
upward shifted by 0.2 eV and 0.3 eV respectively; the energy of the third peak is practically 
unchanged. 

In addition to the above features that originate from the individual bases that enter each 
pair, new features appear in the spectra because the purine and the pyrimidine coexist in the 
Watson-Crick arrangement. This coexistence changes the nature of the frontier orbitals: The 
HOMO is purine-localized, the LUMO is pyrimidine-localized, and the value of the HOMO- 
LUMO gap is smaller than in the isolated bases. Consequently, peaks at lower energies 
emerge in the spectra: they are reported in Table VI, but the corresponding energy range 
is not shown in the figures, because the tiny oscillator strengths make them undiscernible 
from the more intense peaks. We also note that the extremely small intensity makes the 
quantitative assignment of these features less reliable than the others already discussed. In 
the case of GCh, we find excitation energies at 2.44 eV, 3.14 eV and 3.94 eV, with oscillator 
strengths of 0.001, 0.008 and 0.003. These oscillator strengths are two orders of magnitude 
smaller than the peak at 4.29 eV shown in the upper panel of Fig. HI For the ATh pair, 
the lowest energy peak at 4.05 eV has an oscillator strength of 0.048. This transition has a 
small component also in the absorption for light polarized perpendicular to the base plane 
(see bottom right panel in Fig. HJ. Its origin can be likely ascribed to charge transfer states, 
that have received particular attention from both experimental and theoretical viewpoints.— 
However, we warn the reader that the identification of charge transfer transitions is to be 
taken with care at the current level of TDDFT when a local or gradient-corrected exchange- 



20 



HOMO 



LUMO 



GCh 



ATh 




Figure 5. HOMO and LUMO Kohn-Sham wave functions of tlie GCh and ATh pairs, 
correlation functional is used.'^>2Si21 

5. STACKED GC^ AND D(GC) STRUCTURES 

With the results reported in the previous section for the GCh pair, we can understand 
what happens in the optical absorption when a guanine and a cytosine dimerize through 
hydrogen bonding. Here, the GCg system is selected to check which effects arise when the 
same two bases dimerize in a different form, namely by stacking. 

Figure El reports the averaged and direction-resolved spectra of the stacked GCs dimer. 
The comparison between Figs. El and IH aims at disclosing the relative role of hydrogen- 
bonding and TT-stacking in the optical response of nucleobase complexes. We also plot in 
Fig. El the hnear combination of the spectra of the isolated constituent bases. Figure [71 (top) 
illustrates the Kohn-Sham highest occupied and lowest unoccupied orbitals of the system. It 
is important to note that both the HOMO and the LUMO have charge density distributed 
around both the purine and the pyrimidine, at odds with the H-bonded GCh pair (compare 
Fig. El top) where the HOMO is centered at the purine and the LUMO at the pyrimidine. 

As in the case of hydrogen-bonding, the shape of the spectrum it is not strongly altered, 
but important differences are encountered in the oscillator strengths. Looking at the low- 
energy peaks, we see that the excitations due to the cytosine moiety are slightly blue- 
shifted (as for the GCh pair), while the guanine excitations appear unaffected. The resulting 
spectrum below 5eV is made of two equally intense closely spaced peaks at 4.45 eV and 
4.68 eV, whereas in the H-bonded GCh structure we observe two well separated peaks spaced 
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Figure 6. Photo-absorption cross section of the GCs stacked dimer (up) and of the 
d(GC) quartet (down). In each panel, the linear combination of the spectra of the two 
(up) and four (down) constituent bases is also shown. In the right panels the spectra 
resolved in the three spatial directions are shown. 



TT stacked base 



HOMO 



LUMO 



^^^^^^ 



vr stacked pairs '^^3 

Figure 7. Isosurface plots of the HOMO and LUMO Kohn-Sham wave functions of the 
GCs vr-stacked dimer and of the d(GC) 7r-stacked quartet. 



by 0.4 eV and with different intensities. Hypochromicity is also observed in this low-energy 
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Figure 8. Two different conformations obtained for tlie vr-staked GCs pair obtained by 

different structural relaxation schemes. 

range, whereas for the GCh pair hypochromicity occurs only in the high-energy range. 

In the range between 5 and 6eV, we find an enhancement of the weak peak at 5.31 eV 
and a depression of the peak at 5.78 eV, which is present both in the isolated bases (guanine 
and cytosine) and in the H-bonded pair. In the high-energy spectral range, the dominant 
effect is the 39% hypochromicity for the main peak located at 6.15 eV: hypochromicity from 
vr-stacking is thus enhanced with respect to hypochromicity from hydrogen-bonding (30%). 
Above the strongest peak, we find other less intense peaks at 6.47eV, 6.87 eV, and 7.10 eV. 
Instead, in the H-bonded GCh pair, the high-energy range contains a peak at 6.78 eV with 
the same intensity as the peak at 6.17eV, and then only an additional weaker one. 

The calculation of the GCg has been repeated for a geometry obtained by a slightly differ- 
ent relaxation scheme, within the DFT-PW91 framework. ^^'^^ This additional GCg geometry 
was considered just for a gross check of structural deformations induced by computational 
details. The starting configuration was fixed by taking guanine and cytosine from a previous 
DFT study of isolated bases* and putting them parallel to each other at an initial distance of 
3.4 A. The major difference between the two GCs structures (Fig. IH} is the tilt angle, which 
is smaller in the PW91 relaxation. The two spectra are quite similar in both shape and 
intensity: minor differences are observed only in the the first peak, that is shifted upward 
by 0.15 eV in the structure shown in the right panel of Fig. |H1 

The combined effects of stacking and H-bonding on the absorption spectra are simultane- 
ously present in Fig. IHl (lower panel), that corresponds to two stacked GCh pairs with the G's 
on the same strand of the parent polymer.^^ The shape of the spectrum appears very similar 
to that of a single GCh pair (Fig. HI top left), but the intensity of all peaks is reduced. This 
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decrease is a typical qualitative effect of the stacking arrangement, although the molecular 
details fix the quantitative aspects (we described the stacking between guanine and cytosine 
above, whereas we are now focusing on the stacking between two guanines). In fact, in the 
low-energy range we find again very slight shifts for the first two peaks, downward for the 
first and upward for the second, as in Fig. 0] (top left). However, the intensity of the second 
peak is equal to that obtained from the G+C combination, whereas in the single GCh pair 
it was higher. This is due to a depression of the second peak (guanine) induced by stacking. 

In the high-energy range, the hypochromicity due to both the vr-stacking and the H- 
bonding couplings is confirmed. For the main peak located at 6.25 eV we find an intensity 
reduction by 62% with respect to the free monomers. The two kinds of base couplings seem 
to act separately and independently, in the sense that one does not affect the other: such 
separation was already reported concerning the electronic structure.— 

Summarizing this section, we have shown that stacked and H-bonded GC pairs present 
slight differences in the absorption spectra, both in the low-energy and high-energy ranges. 
Hypochromicity has been found in both configurations and is largest for the stacked pairs. 
Hypochromicity is very useful because this intensity change can be used to follow the melting 
of the secondary structure of nucleic acids when varying the temperature or environmental 
parameters. 

6. SUMMARY AND PERSPECTIVES 

In this work we presented a complete study of the optical absorption spectra of the five 
isolated gas-phase nucleobases and their assemblies: simple Watson-Crick pairs, simple vr- 
stacks of two bases and more complex vr-stacks of Watson-Crick GC pairs. These calculations 
of isolated simple nucleobase-assemblies are a remarkable playground to prove the reliability 
of TDDFT for DNA-based materials. Remarkably, for the first time the optical properties 
were computed by an ab-initio method for a helical conformation of two stacked GC pairs, 
where H-bonding and vr-stacking effects are active simultaneously and can be distinguished. 

The results can be summarized as follows. For the isolated bases we get spectra in good 
agreement with previous theoretical works and (qualitatively) with experiments. We repro- 
duce the proper ordering of the vrvr* excitations, namely the excitation energy increases in 
going from C to G to A. Moreover, the LUMO state has always a vr-like character whereas 
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the HOMO is vr-like for the purines and cr-hke for the pyrimidines. As concerns the base 
assembhes (Watson-Crick H-bond pairs and stacked configuration), we obtain that the shape 
of the spectrum is not much altered by the vr-stacking or H-bond interactions. However, we 
always get hypochromicity in the high energy range of the spectrum. The hypochromicity 
induced by vr-stacking is larger than that induced by H-bonding. For light polarized perpen- 
dicular to the bases we get a blue-shift of the spectra compared to the spectra of the isolated 
bases. In the stacked case, the HOMO and LUMO states are distributed both on the purine 
and pyrimidine bases, whereas in the H-bonded configuration the HOMO is in the purine 
and the LUMO in the pyrimidine (charge transfer- type excitation). When combining both 
H-bonding and vr-stacking, the two effects add independently, and the hypochromicity in 
the UV is enhanced. 

At this stage it is relevant to note that all our calculations are in the gas phase, which 
means that solvent and environment effects are not taken into account. This hinders a 
direct comparison with experimental data, as most of the available results correspond to 
nucleobases assemblies in solution. However, several experimental and theoretical studies 
(including simplified models of the solvent) showed that the vrvr* transitions are only slightly 
affected by the presence of the solven1)^>^2^: precisely, they are insensitive to the polarity 
of the solvent, but may undergo red-shifts. On the contrary, the "dark" nvr* transitions 
turn out to be significantly sensible to the polarity of the solvent . ^^i^^i^'^'i^^i^^i'^^ As shown 
in Figs. 121 and the absorption spectra of the DNA-bases is extremely anisotropic, due to 
their quasi-planar structure. Consequently, nearly all the oscillator strength is concentrated 
on the vrvr* transitions (orders of magnitude more intense than the nvr* transitions). Those 
transitions are excited only for light polarized in the base plane. Therefore, these results 
can be used to discriminate (albeit qualitatively) which features are mostly limited to the 
gas phase (those polarized out-of-plane) and which ones can be considered intrinsic to the 
DNA-complex.iSii Furthermore, the role of charge-transfer states needs of be analyzed more 
in detail in order to understand their impact in the excited state dynamics of DNA-based 
complexes (note that those states are very likely to be dark, or with very low oscillator 
strength for light-induced electronic excitations). 

Besides optical absorption, there is another optical technique that is widely used for the 
characterization of chiral biomolecules - circular dichroism. The computational method used 
in the present work allows for a straightforward calculation of the rotatory power or circular 
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dichroism spectra keeping the simplicity of the time-propagation scheme, i.e, without the 
need for empty states (see Refj4^for the details). The implementation and computation of 
circular dichroism spectra by TDDFT is in progress and will be the topic of a self-standing 
investigation,^^ that should allow the identification of helical fingerprint in the optical char- 
acteristics, and more direct interpretation of standard post-synthesis experimental data. 

Acknowledgements 

The authors thank X. Lopez and E. Sansebastian for helpful discussion. Funding was 
provided by the EC through grant IST-200 1-38951, by the EC 6th Framework Network of 
Excellence NANOQUANTA (NMP4-CT-2004-500198), by Spanish MFC and by the INFM 
through the Parallel Computing Committee. A. Rubio acknowledges the Humboldt Foun- 
dation under the Bessel research award (2005). We thank computer time provided by the 
SGI of the UPV and DIPC as well as the Barcelona supercomputer center. 

References and Notes 

^ Daniels, M. In Photochemistry and Photobiology of Nucleic Acids: Wang S. Y. Ed; 
Academic Press: New York, 1976; p. 23. 

2 Hernandez, C. E. C; Cohen, B.; Hare, P. M.; Kohler, B. Chem. Rev. 2004, 104, 1977. 

^ Ruzsicska, B. P.; Lemaire, D. G. E. In CRC Handbook of Photochemistry and Photo- 
biology; HorspoU W. M.; Lemaire D. G. E. Eds; CRC Press: Boca Raton, Fl, 1995; 
p. 1289 

^ Porath, D.; Cuniberti, G.; Felice, R. D. In Topics in Current Chemistry; Schuster G. 
Ed.; Springer Verlag, Heideleberg, 2004; p. 183. 

^ Di Ventra, M.; Zwolak, M. In Encycolpedia of Nanoscience and Nanotechnology and 
Handbook of Nanostructured Biomaterials and their application in Nanotechnology; 
Nalwa H. S. Ed.; APS, 2004; p 475. 

6 Endres, R. G.; Cox, D.; Singh, R. Rev. Mod. Phys. 2004, 76, 195. 



26 



^ Coehn, H.; Nogues, C; Naaman, R.; Porath, D. Proc. Natl. Acad. Sci. U.S.A. 2005, 
102, 11589. 

^ Di Felice, R.; Calzolari, A.; Molinari, E.; Garbesi, A. Phys. Rev. B 2002, 65, 045104. 

^ Di Felice, R.; A. Calzolari. In Modem Methods for theoretical physical chemestry of 
biopolymers; Starikov, E. B.; Tanaka, S.; Lewis J. R Eds; Elsevier, 2005. 

Kraemer, K. H. Proc. Natl. Acad. Sci. U.S.A. 1997, 94, 11. 

^ Odom, D. T.; Barton, J. K. Biochemistry 2001, 40, 8727. 

2 Hernandez, C. E. C; Kohler, B.; Cohen, B. Nature 2005, 436, 1141. 

3 Watson, J. H.; Crick, F. H. Nature 1953, 171, 737. 

^ Hall, D. B.; Holmlin, R. E.; Barton, J. K. Nature 1996, 3841, 731. 
^ Runge, E.; Gross, E. K. U. Phys. Rev. Lett. 1984, 52, 997. 

^ Gross, E. K. U.; Dobson, J. F.; Petersilka, M. Density functional theory II, Topics in 
Current Chemistry; Nalewajski Ed; Springer Verlag: Berlin, 1996. 

^ Marques, M. A. L.; Gross, E. K. U. In A primer in Density Functional Theory; Fiolhais 
C; Nogueira F.; Marques M. A. L Marques Eds; Springer: Berlin, 2003: p.l44. 

s Onida, G.; Reining, L.; Rubio, A. Rev. Mod. Phys. 2002, 74, 601. 

^ Furche, F.; Burke, K. Ann. Rev. Comp. Chem. 2005, 1, 19. 

'° Marques, M. A. L.; Lopez, X.; D.Varsano.; Castro, A.; Rubio, A. Phys. Rev. Lett. 
2003, 90, 258101. 

Lopez, X.; Marques, M. A. L.; Castro, A.; Rubio, A. J. Am. Chem. Soc. 2005, 127, 
12329. 

'2 Sundholm, D. Phys. Chem. Chem. Phys. 2000, 2, 2275. 
Sundholm, D. Phys. Chem. Chem. Phys. 2003, 5, 4265. 

Neiss, C; Saalfrank, R; Parac, M.; Grimme, S. J. Phys. Chem. A 2003, 107, 140. 

27 



25 Tsolakidis, A.; Kaxiras, E. J. Phys. Chem. A 2005, 109, 2373. 

2^ Shukla, M. K.; Leszczynski, J. J. Comput. Chem. 20 04, 25, 768. 

2^ von Lilienfeld, O. A.; Tavernelli, I.; Rothlisberger, U.; Sebastiani, D. Phys. Rev. Lett. 
2004, 93, 153004. 



28 



Lorentzon, J.; Fulscher, M. P.; Roos, B. J. Am. Chem. Soc. 1995, 117, 9265. 



2^ Fulscher, M. R; Serrano-Andres, L.; Roos, B. J. Am. Chem. Soc. 1997, 119, 6168. 
3° Fulscher, M. P.; Roos, B. J. Am. Chem. Soc. 1995, 117, 2089. 
3^ Shukla, M. K.; Leszczynski, J. J. Phys. Chem. A 2002, 106, 1011. 

32 Shukla, M. K.; Leszczynski, J. J. Phys. Chem. A 20 02, 106, 4709. 

33 Sobolewski, A.; Domcke, W. Phys. Chem. Chem. Phys. 2004, 6, 2763. 

34 Wesoloswki, T. J. Am. Chem. Soc. 2004, 126, 11444. 

35 Jean, J. M.; Hall, K. B. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 37. 

36 Jean, J. M.; Hall, K. B. Biochemistry 2002, 4I, 13152. 

3^ Mennucci, B.; Tonniolo, A.; Tommasi, J. J. Phys. Chem. A 2001, 105, 4749. 

38 Mishra, S. K.; Shukla, M. K.; Mishra, P. C. Spectrochim. Acta Part A 2000, 56, 1355. 

39 Shukla, M. K.; Mishra, S. K.; Kumar, A.; Mishra, P. C. J. Comput. Chem. 2000, 21, 
826. 

4° Shukla, M. K.; Leszczynski, J. J. Phys. Chem. A 2002, 106, 11338. 

Marques, M. A. L.; Castro, A.; Bertsch, G. F.; Rubio, A. Computer Phys. Commun. 
2003, 151, 60; available at|http:// www.tddft.org/programs/octopusl 



42 Castro, A.; Marques, M. A. L.; Alonso, J.; Rubio, A. J. of Comput. and Theor. 
Nanoscience 2004, 1, 230. 

43 Yabana, K.; Bertsch, G. F. Phys. Rev. B 1996, 5^, 4484. 



28 



^ Yabana, K.; Bertsch, G. F. Int. J. Quantum Chem. 1999, 75, 55. 

^5 Yabana, K.; Bertsch, G. F. Phys. Rev. A 1999, 60, 3809. 

Rubio, A.; Alonso, J.; Blase, X.; .; Balbas, L.; Louie, S. Phys. Rev. Lett. 1996, 77, 247. 

Rubio, A.; Alonso, J.; Blase, X.; Louie, S. International Journal of Modern Physics B 
1997, 11, 2727. 

Casida, E. In Recent Advances in Density Junctionals Methods Part 1; Chong D. P Ed.; 
World Scientific: Singapore, 1995; p. 115. 

Perdew, J.; Zunger, A. Phys. Rev. B 1981, 23, 5048. 

Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865. 

^1 Marques, M. A. L.; Castro, A.; Rubio, A. J. Chem Phys. 2001, 115, 3006. 

In spite of the equivalence encountered in the photo- absorption cross-section of gas- 
phase guanine by using these two functionals, which is sufficient for our purposes to jus- 
tify the LDA option, one has to be extremely careful in general in the choice of Vxc- For 
instance, it was shown that in processes that involve charge transfer mechanisms, like 
fluorescence in base pairs, the correct asymptotic behavior of the exchange-correlation 
functional may affect the accuracy of the outcome. See, e.g., S. Grimme and M. Parac, 
ChemPhysChem 2003, ^, 292 . 

53 TrouUier, N.; Martins, J. L. Phys. Rev. B 1991, 43, 1993. 

5^ Arnott, S.; Hukins, D. W. L. Biochem. Biophys. Res. Commun. 1972, 47, 1504. 

55 Gao, Y. G.; Robinson, H.; Wang, A. H. Eur. J. Biochem. 1999, 261, 413. 

56 Brunger, A. T.; Karplus, M. Proteins 1988, 4, 148. 

57 Brooks, B. R.; others. J. Comput. Chem. 1983, 2, 148. 

58 MacKerrel, A. D.; others. J. Phys. Chem. B 1998, 102, 3586. 

59 Perdew, J. P.; Wang, Y. Phys. Rev. B 1992, 45, 13244. 
Zhang, H.; Di Felice, R. unpublished. 

29 



Calzolari, A.; Di Felice, R.; Molinari, E.; Garbesi, A. J. Phys. Chem. B 2004, 108, 
2509. 

^2 Calzolari, A.; Di Felice, R.; Molinari, E.; Garbesi, A. J. Phys. Chem. B 2004, 108, 
13058. 

Shukla, M. K.; Mishra, S. K. Chem. Phys. 1999, 2^0, 319. 

^ Merchan, M.; Serrano- Andres, L. J. Am. Chem. Soc. 2003, 125, 8108. 

Ismail, N.; Blancafort, L.; Olivucci, M.; Kohler, B.; Robb, M. A. J. Am. Chem. Soc. 
2002, 87, 6818. 

^6 Clark, L. B.; Tinoco, I. J. Am. Chem. Soc. 1965, 87, 11. 

6^ Miles, D. W.; Robins, R. K.; Eyring, K. Proc. Natl. Acad. Sci. U.S.A. 1967, 57, 1139. 

Sprecher, C. A.; Johnson Jr, W. C. Biopolymers 1977, 16, 2243. 

Voelter, W.; Records, R.; Bunnenberg, E.; Djerassi, C. J. Am. Chem. Soc. 1968, 90, 
6163. 

™ Fucaloro, A. F.; Forster, L. S. J. Am. Chem. Soc. 1971, 93, 6443. 

Eaton, W. A.; Lewis, T. R J. Chem. Phys. 1970, 53, 2164. 

^2 Stewart, R. F.; Davidson, N. J. Chem. Phys. 1963, 39, 255. 

'^^ Stewart, R. F.; Davidson, N. Biopolymers 1964, 1, 465. 

Kim, N. J.; Jeong, G.; Kim, Y. S.; Sung, J.; Kim, S. K.; Park, Y. J. Chem. Phys. 
2000, 113, 10051. 

^5 Sobolewski, A.; Domcke, W. Eur. Phys. J. D 2002, 20, 369. 

^6 Holmen, A.; Broo, A.; Albinsson, B.; Norden, B. J. Am. Chem. Soc. 1997, 119, 12240. 
Broo, A. J. Phys. Chem. A 1998, 102, 526. 

Callis, R R.; Rosa, E. J.; Simpson, W. T. J. Am. Chem. Soc. 1964, 86, 2292. 
^9 Clark, L. B. J. Phys. Chem. 1995, 99, 4466. 

30 



Brunner, W. C; Maestre, M. F. Biopolymers 1975, I4, 555. 
Hernandez, C. E. C; Marai, C; Kohler, B. unpublished 2004. 
Nir, E.; others. Chem. Phys. Lett. 2002, 59, 355. 

Mons, M.; Dimicoli, I.; F.Piuzzi.; Tardivel, B.; Elhanine, M. J. Phys. Chem. A 2002, 
106, 5088. 

Jortner, J; Bixon, M.; Langenbacher, T; Michel-Beyerle, M. E. Proc. Natl. Acad. Sci. 
U.S.A. 1998, 95, 12759. 

Dreuw, A.; Weisman, J. L.; Head-Gordon, M. J. Chem. Phys. 2003, 119, 2943. 
Dreuw, A.; Weisman, J. L.; Head-Gordon, M. J. Am. Chem. Soc. 2004, 126, 4007. 
Gritsenko, O; Baerends E. J. J. Chem. Phys. 2004, 121, 655. 
Brealey, G.; Kasha, M. J. Am. Chem. Soc. 1955, 77, 4462. 
Calhs, P. R. Annu. Rev. Phys. Chem. 1983, 37, 691. 

With this reasoning, we do not claim that the values of the tttt* transitions are to be 
considered exact also for the solution phase, but they are not expected to undergo 
major modifications in different solutions. Most likely, the effect of the environment on 
the TTTT* transitions can be expected to be a systematic shift toward low energies, that 
should be taken into account for a more strict comparison in Tables I- VI. However, 
such a shift cannot be computed at present at the same level of theory. 

Varsano, D.; Felice, R. D.; Marques, M. A. L.; Rubio, A. unpublished. 



31 



